Signal processing method and base station

ABSTRACT

Embodiments of the present invention provide a signal processing method and a base station. By using respective sparse characteristics of a user signal and an interference signal, the user signal and the interference signal are iteratively recovered, that is, the interference signal is first fastened, and the user signal is recovered; and then the user signal is fastened, and the interference signal is recovered. In this way, an iteration cycle is repeated until the recovered user signal and the recovered interference signal meet a preset condition, thereby effectively improving interference cancellation performance, achieving relatively good interference suppression.

CROSS-REFERENCE TO RELATED APPLICATIONS

This application is a continuation of International Application No. PCT/CN2014/072571, filed on Feb. 26, 2014, which is hereby incorporated by reference in its entirety.

TECHNICAL FIELD

Embodiments of the present invention relate to the field of communications technologies, and in particular, to an information processing method and a base station.

BACKGROUND

Compressive sensing is an emerging signal processing technology in recent years, and a core idea of the compressive sensing is to simultaneously sample and compress data. A non-adaptive linear projection (a measurement value) of a signal is first collected, and then the signal is recovered by using the measurement value and according to a corresponding reconstruction algorithm. The compressive sensing has basic requirements in two aspects: sparseness of a signal and a non-correlation between an observation group (an observation matrix) and a transform group (a transform matrix). Any signal in nature has particular representation space, so that the signal has sparseness in this space.

A capacity of a wireless communications system is limited by interference. For example, in a mobile cellular system, there is not only intra-cell multiuser interference, but also inter-cell interference. Therefore, interference control and cancellation is a key technology for ensuring normal operation of a wireless network. An existing wireless interference control technology is mainly divided into three types as follows:

An interference avoidance technology: for example, Frequency Division Multiple Access (Frequency Division Multiple Access, FDMA for short), Time Division Multiple Access (Time Division Multiple Access, TDMA for short), and Orthogonal Frequency Division Multiple Access (Orthogonal Frequency Division Multiple Access, OFDMA for short), and these multiple access technologies are essentially to send, in mutually-orthogonal signal space, signals of different users in neighboring cells, thereby avoiding mutual interference between sent signals.

An interference averaging technology: in Code Division Multiple Access (Code Division Multiple Access, CDMA for short), signals sent by users are all extended by different pseudo-random codes to entire signal space, thereby implementing interference averaging and interference suppression.

A cooperative interference suppression technology: some new interference suppression technologies based on signal processing and wireless node cooperation have emerged recently. For example, interference alignment implements interference suppression by centralizing interference in same sub-space. A cooperative multiple-input multiple-output technology (Cooperative Multiple-Input Multiple-Output, Co-MIMO for short) is used to perform interference cancellation and improve a transmission rate by sharing channel state information and transmitting data between transmit ends that interfere with each other.

However, the foregoing wireless interference control technologies have respective disadvantages. For example, the interference avoidance technology cannot fully use signal latitude, and therefore has relatively low performance; the interference averaging technology has relatively poor performance when a signal-to-noise ratio is high; for the cooperative interference suppression technology, an advantage of the technology cannot be made used of because current wireless communications technologies have a lot of problems that are not resolved. Therefore, none of the foregoing technical solutions can implement optimal interference suppression.

SUMMARY

Embodiments of the present invention provide a signal processing method and a base station, to implement interference suppression, further improving interference cancellation performance, and a signal-to-noise ratio and a diversity gain that are of a receive end.

According to a first aspect, an embodiment of the present invention provides a signal processing method, including:

receiving, by a base station, a signal y, where the signal y includes a user signal X_(S) and an interference signal G_(I), where: X_(S) is a vector including S user signals, and X_(S)=[x₁ x₂ x₃ . . . x_(s)], where x₁ represents a first user signal, x₂ represents a second user signal, x₃ represents a third user signal, and so on, and x_(s) represents an s^(th) user signal; and G_(I) is a vector including I interference signals, and G_(I)=[g₁ g₂ g₃ . . . g_(I)], where g₁ represents a first interference signal, g₂ represents a second interference signal, g₃ represents a third interference signal, and so on, and g_(I) represents an I^(th) interference signal;

determining, by the base station, an initial value G_(I) ⁰ of the interference signal G_(I), where G_(I) ⁰=[g₁ ⁰ g₂ ⁰ g₃ ⁰ . . . g_(I) ⁰];

determining, by the base station, a value X_(S) ¹ of the user signal X_(S) according to the initial value G_(I) ⁰ and the signal y, where X_(S) ¹=[x₁ ¹ x₂ ¹ x₃ ¹ . . . x_(s) ¹];

determining, by the base station, that the value X_(S) ¹ is an initial value X_(S) ⁰ of the user signal X_(S), where X_(S) ⁰=[x₁ ¹ x₂ ¹ x₃ ¹ . . . x_(s) ¹];

determining, by the base station, a value G_(I) ¹ of the interference signal G_(I) according to the initial value X_(S) ⁰ and the signal y, where G_(I) ¹=[g₁ ¹ g₂ ¹ g₃ ¹ . . . g_(I) ¹];

determining, by the base station, that the value G_(I) ¹ is the initial value G_(I) ⁰ of the interference signal G_(I), where G_(I) ⁰=[g₁ ¹ g₂ ¹ g₃ ¹ . . . g_(I) ¹];

performing, by using the foregoing method, iterative processing on the determined value X_(S) ¹ of the user signal X_(S) and the determined value G_(I) ¹ of the interference signal G_(I), until the value X_(S) ¹ of the user signal X_(S) and the value G_(I) ¹ of the interference signal G_(I) that are obtained through an m^(th) iteration meet a preset condition, determining that an approximate value of the user signal X_(S) is X_(S) ¹, and determining that an approximate value of the interference signal G_(I) is G_(I) ¹.

In a first possible implementation manner of the first aspect, the determining, by the base station, an initial value G_(I) ⁰ of the interference signal G_(I) includes:

solving, by the base station, a formula (1) to obtain Z, where Z includes an optimal solution G_(I) ^(opt) of the interference signal G_(I):

$\begin{matrix} {{\min\limits_{z}{{{\Phi \; Z} - y}}_{2}},{{s.t.{Z}_{0}} = {S + I}}} & (1) \end{matrix}$

where ∥ ∥₂ denotes the l₂-norm; ∥ ∥₀ denotes the l₀-norm; Φ=[Φ_(S) Φ_(I)], and

${Z = \begin{bmatrix} X_{S} \\ G_{I} \end{bmatrix}};$

Φ_(S)=H or Φ_(S)=H_(S)f_(S), where H_(s) is a frequency-domain response of a channel passed through by the user signal, and f_(s) is a precoding vector used by user equipment UE; and Φ_(I)=E_(I) or Φ_(I)=H_(I)f_(I), where E_(I) is a unit vector, H_(I) is a frequency-domain response of a channel passed through by the interference signal, and f_(I) is a precoding vector used by interference user equipment UE; and

determining, by the base station, that the optimal solution G_(I) ^(opt) is the initial value G_(I) ⁰ of the interference signal G_(I).

In a second possible implementation manner of the first aspect, the determining, by the base station, a value X_(S) ¹ of the user signal X_(S) according to the initial value G_(I) ⁰ and the signal y includes:

solving, by the base station, a formula (2) to obtain an optimal solution Z_(S), and determining that Z_(S) is the value X_(S) ¹ of the user signal X_(S):

$\begin{matrix} {{\min\limits_{z_{S}}{{{\Phi_{S}\; Z_{S}} - y_{S}}}_{2}},{{s.t.{Z_{S}}_{0}} = S}} & (2) \end{matrix}$

where y_(S)=y−Φ_(I)G_(I) ⁰; and

correspondingly, the determining, by the base station, a value G_(I) ¹ of the interference signal G_(I) according to the initial value X_(S) ⁰ and the signal y includes:

solving, by the base station, a formula (3) to obtain an optimal solution Z_(I), and determining that Z_(I) is the value G_(I) ¹ of the interference signal G_(I):

$\begin{matrix} {{\min\limits_{z_{I}}{{{\Phi_{I}\; Z_{I}} - y_{I}}}_{2}},{{s.t.{Z_{I}}_{0}} = I}} & (3) \end{matrix}$

where y_(I)=y−Φ_(S)X_(S) ⁰.

With reference to the first or the second possible implementation manner of the first aspect, in a third possible implementation manner of the first aspect, the base station solves a formula by using a compressive sensing signal recovery algorithm to obtain an optimal solution.

With reference to the first aspect, the first possible implementation manner of the first aspect, the second possible implementation manner of the first aspect, or the third possible implementation manner of the first aspect, in a fourth possible implementation manner of the first aspect, the performing, by the base station by using the foregoing method, iterative processing on the determined value X_(S) ¹ of the user signal X_(S) and the determined value G_(I) ¹ of the interference signal G_(I), until the value X_(S) ¹ of the user signal X_(S) and the value G_(I) ¹ of the interference signal G_(I) that are obtained through an m^(th) iteration meet a preset condition includes:

performing, by the base station by using the foregoing method, iterative processing on the determined value X_(S) ¹ of the user signal X_(S) and the determined value G_(I) ¹ of the interference signal G_(I), until a difference between an estimated error ε^(m-1) obtained through calculation of a formula (4) by using the value X_(S) ¹ of the user signal X_(S) and the value G_(I) ¹ of the interference signal G_(I) that are obtained through an (m-1)^(th) iteration and an estimated error ε^(m) obtained through calculation of the formula (4) by using the value X_(S) ¹ of the user signal X_(S) and the value G_(I) ¹ of the interference signal G_(I) that are obtained through the m^(th) iteration is greater than a preset threshold ε:

ε^(i)=∥Φ_(S) X _(S) ¹+Φ_(I) G _(I) ¹ −y∥₂   (4)

where i is a natural number, and ε^(i) represents an estimated error obtained through an i^(th) iteration;

or

collecting, by the base station, statistics on a quantity of iteration times, until the quantity of iteration times is greater than a preset threshold m.

According to a second aspect, an embodiment of the present invention provides a base station, including:

a receiving module, configured to receive a signal y, where the signal y includes a user signal X_(S) and an interference signal G_(I), where: X_(S) is a vector including S user signals, and X_(S)=[x₁ x₂ x₃ . . . x_(s)], where x₁ represents a first user signal, x₂ represents a second user signal, x₃ represents a third user signal, and so on, and x_(s) represents an s^(th) user signal; and G_(I) is a vector including I interference signals, and G_(I)=[g₁ g₂ g₃ . . . g_(I)], where g₁ represents a first interference signal, g₂ represents a second interference signal, g₃ represents a third interference signal, and so on, and g_(I) represents an I^(th) interference signal; and

a determining module, configured to determine an initial value G_(I) ⁰ of the interference signal G_(I), where G_(I) ⁰=[g₁ ⁰ g₂ ⁰ g₃ ⁰ . . . g_(I) ⁰]; determine a value X_(S) ¹ of the user signal X_(S) according to the initial value G_(I) ⁰ and the signal y, where X_(S) ¹=[x₁ ¹ x₂ ¹ x₃ ¹ . . . x_(s) ¹]; determine that the value X_(S) ¹ is an initial value X_(S) ⁰ of the user signal X_(S), where X_(S) ⁰=[x₁ ¹ x₂ ¹ x₃ ¹ . . . x_(s) ¹]; determine a value G_(I) ¹ of the interference signal G_(I) according to the initial value X_(S) ⁰ and the signal y, where G_(I) ¹=[g₁ ¹ g₂ ¹ g₃ ¹ . . . g_(I) ¹]; determine that the value G_(I) ¹ is the initial value G_(I) ⁰ of the interference signal G_(I), where G_(I) ⁰=[g_(I) ¹ g_(I) ² g_(I) ³ . . . g_(I) ¹]; and perform iterative processing on the value X_(S) ¹ of the user signal X_(S) and the value G_(I) ¹ of the interference signal G_(I), until the value X_(S) ¹ of the user signal X_(S) and the value G_(I) ¹ of the interference signal G_(I) that are obtained through an m^(th) iteration meet a preset condition, determine that an approximate value of the user signal X_(S) is X_(S) ¹, and determine that an approximate value of the interference signal G_(I) is G_(I) ¹.

In a first possible implementation manner of the second aspect, the determining module is specifically configured to:

solve a formula (1) to obtain Z, where Z includes an optimal solution G_(I) ^(opt) of the interference signal G_(I):

$\begin{matrix} {{\min\limits_{z}{{{\Phi \; Z} - y}}_{2}},{{s.t.{Z}_{0}} = {S + I}}} & (1) \end{matrix}$

where ∥ ∥₂ denotes the l₂-norm; ∥ ∥₀ denotes the l₀-norm; Φ=[Φ_(S) Φ_(I)], and

${Z = \begin{bmatrix} X_{S} \\ G_{I} \end{bmatrix}};$

Φ_(S)=H or Φ_(S)=H_(S)f_(S), where H_(s) is a frequency-domain response of a channel passed through by the user signal, and f_(s) is a precoding vector used by user equipment UE; and Φ_(I)=E_(I) or Φ_(I)=H_(I)f_(I), where E_(I) is a unit vector, H_(I) is a frequency-domain response of a channel passed through by the interference signal, and f_(I) is a precoding vector used by interference user equipment UE; and determine that the optimal solution G_(I) ^(opt) is the initial value G_(I) ⁰ of the interference signal G_(I).

In a second possible implementation manner of the second aspect, the determining module is specifically configured to:

solve a formula (2) to obtain an optimal solution Z_(S), and determine that Z_(S) is the value X_(S) ¹ of the user signal X_(S):

$\begin{matrix} {{\min\limits_{z_{S}}{{{\Phi_{S}\; Z_{S}} - y_{S}}}_{2}},{{s.t.{Z_{S}}_{0}} = S}} & (2) \end{matrix}$

where y_(S)=y−Φ_(I)G_(I) ⁰; and solve a formula (3) to obtain an optimal solution Z_(I), and determine that Z_(I) is the value G_(I) ¹ of the interference signal G_(I):

$\begin{matrix} {{\min\limits_{z_{I}}{{{\Phi_{I}\; Z_{I}} - y_{I}}}_{2}},{{s.t.{Z_{I}}_{0}} = I}} & (3) \end{matrix}$

where y_(I)=y−Φ_(S)X_(S) ⁰.

With reference to the first or the second possible implementation manner of the second aspect, in a third possible implementation manner of the second aspect, the base station further includes a compressive sensing module, configured to solve, by using a compressive sensing signal recovery algorithm, a formula to obtain an optimal solution.

With reference to the second aspect, the first possible implementation manner of the second aspect, the second possible implementation manner of the second aspect, or the third possible implementation manner of the second aspect, in a fourth possible implementation manner of the second aspect, the determining module is specifically configured to:

perform iterative processing on the value X_(S) ¹ of the user signal X_(S) and the value G_(I) ¹ of the interference signal G_(I), until a difference between an estimated error ε^(m-1) obtained through calculation of a formula (4) by using the value X_(S) ¹ of the user signal X_(S) and the value G_(I) ¹ of the interference signal G_(I) that are obtained through an (m-1)^(th) iteration and an estimated error ε^(m) obtained through calculation of the formula (4) by using the value X_(S) ¹ of the user signal X_(S) and the value G_(I) ¹ of the interference signal G_(I) that are obtained through the m^(th) iteration is greater than a preset threshold ε:

ε^(i)=∥Φ_(S) X _(S) ¹+Φ_(I) G _(I) ¹ −y∥₂   (4)

where i is a natural number, and ε^(i) represents an estimated error obtained through an i^(th) iteration; or collect statistics on a quantity of iteration times, until the quantity of iteration times is greater than a preset threshold m.

According to the signal processing method and the base station provided in the embodiments of the present invention, by using respective sparse characteristics of a user signal and an interference signal, the user signal and the interference signal are iteratively recovered by using a compressive sensing signal recovery algorithm, that is, the interference signal is first fastened, and the user signal is recovered by using the compressive sensing signal recovery algorithm; and then the user signal is fastened, and the interference signal is recovered by using the compressive sensing signal recovery algorithm. In this way, an iteration cycle is repeated until the recovered user signal and the recovered interference signal meet a preset condition, thereby effectively improving interference cancellation performance, achieving relatively good interference suppression.

BRIEF DESCRIPTION OF DRAWINGS

To describe the technical solutions in the embodiments of the present invention or in the prior art more clearly, the following briefly introduces the accompanying drawings required for describing the embodiments or the prior art. Apparently, the accompanying drawings in the following description show some embodiments of the present invention, and persons of ordinary skill in the art may still derive other drawings from these accompanying drawings without creative efforts.

FIG. 1 is a schematic flowchart of Embodiment 1 of a signal processing method according to the present invention;

FIG. 2 is a schematic diagram of an application scenario 1 according to the present invention;

FIG. 3 is a schematic diagram of precoding processing in the application scenario 1 according to the present invention;

FIG. 4 is a schematic diagram of an application scenario 2 according to the present invention;

FIG. 5 is a schematic flowchart of Embodiment 2 of a signal processing method according to the present invention;

FIG. 6 is a schematic flowchart of Embodiment 3 of a signal processing method according to the present invention; and

FIG. 7 is a schematic structural diagram of an embodiment of a base station according to the present invention.

DESCRIPTION OF EMBODIMENTS

To make the objectives, technical solutions, and advantages of the embodiments of the present invention clearer, the following clearly and completely describes the technical solutions in the embodiments of the present invention with reference to the accompanying drawings in the embodiments of the present invention. Apparently, the described embodiments are some but not all of the embodiments of the present invention. All other embodiments obtained by persons of ordinary skill in the art based on the embodiments of the present invention without creative efforts shall fall within the protection scope of the present invention.

Embodiment 1

FIG. 1 is a schematic flowchart of Embodiment 1 of a signal processing method according to the present invention, and the signal processing method specifically includes the following steps:

S101: A base station receives a signal y, where the signal y includes a user signal X_(S) and an interference signal G_(I), where: X_(S) is a vector including S user signals, and X_(S)=[x₁ x₂ x₃ . . . x_(s)], where x₁ represents a first user signal, x₂ represents a second user signal, x₃ represents a third user signal, and so on, and x_(s) represents an s^(th) user signal; and G_(I) is a vector including I interference signals, and G_(I)=[g₁ g₂ g₃ . . . g_(I)], where g₁ represents a first interference signal, g₂ represents a second interference signal, g₃ represents a third interference signal, and so on, and g_(I) represents an I^(th) interference signal.

S102: The base station determines an initial value G_(I) ⁰ of the interference signal G_(I), where G_(I) ⁰=[g₁ ⁰ g₂ ⁰ g₃ ⁰ . . . g_(I) ⁰].

Specifically, the base station may determine the initial value G_(I) ⁰ of the interference signal G_(I) by using the following method:

solving, by the base station, a formula (1) to obtain Z, where Z includes an optimal solution G_(I) ^(opt) of the interference signal G_(I):

$\begin{matrix} {{\min\limits_{z}{{{\Phi \; Z} - y}}_{2}},{{s.t.{Z}_{0}} = {S + I}}} & (1) \end{matrix}$

where ∥ ∥₂ denotes the l₂-norm; ∥ ∥₀ denotes the l₀-norm; Φ=[Φ_(S) Φ_(I)], and

${Z = \begin{bmatrix} X_{S} \\ G_{I} \end{bmatrix}};$

Φ_(S)=H or Φ_(S)=H_(S)f_(S), where H_(s) is a frequency-domain response of a channel passed through by the user signal, and f_(s) is a precoding vector used by user equipment UE; and Φ_(I)=E_(I) or Φ_(I)=H_(I)f_(I), where E_(I) is a unit vector, H_(I) is a frequency-domain response of a channel passed through by the interference signal, and f_(I) is a precoding vector used by interference user equipment UE; and determining that the optimal solution G_(I) ^(opt) is the initial value G_(I) ⁰ of the interference signal G_(I).

For example, the base station may solve, by using a compressive sensing signal recovery algorithm, the foregoing formula (1) to obtain the optimal solution, for example, a BP algorithm, an MP algorithm, an OMP algorithm, an ROMP algorithm, a CoSaMP algorithm, and an SP algorithm.

S103: The base station determines a value X_(S) ¹ of the user signal X_(S) according to the initial value G_(I) ⁰and the signal y, where X_(S) ¹=[x₁ ¹ x₂ ¹ x₃ ¹ . . . x_(s) ¹].

Specifically, the base station solves a formula (2) to obtain an optimal solution Z_(S), and determines that Z_(S) is the value X_(S) ¹ of the user signal X_(S):

$\begin{matrix} {{\min\limits_{z_{S}}{{{\Phi_{S}\; Z_{S}} - y_{S}}}_{2}},{{s.t.{Z_{S}}_{0}} = S}} & (2) \end{matrix}$

where y_(S)=y−Φ_(I)G_(I) ⁰.

For example, the base station may solve, by using the compressive sensing signal recovery algorithm, the foregoing formula (2) to obtain the optimal solution, for example, a BP algorithm, an MP algorithm, an OMP algorithm, an ROMP algorithm, a CoSaMP algorithm, and an SP algorithm.

S104: The base station determines that the value X_(S) ¹ is an initial value X_(S) ⁰ of the user signal X_(S), where X_(S) ⁰=[x₁ ¹ x₂ ¹ x₃ ¹ . . . x_(s) ¹].

S105: The base station determines a value G_(I) ¹ of the interference signal G_(I) according to the initial value X_(S) ⁰ and the signal y, where G_(I) ¹=[g₁ ¹ g₂ ¹ g₃ ¹ . . . g_(I) ¹].

Specifically, the base station solves a formula (3) to obtain an optimal solution Z_(I), and determines that Z_(I) is the value G_(I) ¹ of the interference signal G_(I):

$\begin{matrix} {{\min\limits_{z_{I}}{{{\Phi_{I}\; Z_{I}} - y_{I}}}_{2}},{{s.t.{Z_{I}}_{0}} = I}} & (3) \end{matrix}$

where y_(I)=y−Φ_(S)X_(S) ⁰.

For example, the base station may solve, by using the compressive sensing signal recovery algorithm, the foregoing formula (2) to obtain the optimal solution, for example, a BP algorithm, an MP algorithm, an OMP algorithm, an ROMP algorithm, a CoSaMP algorithm, and an SP algorithm.

S106: The base station determines that the value G_(I) ¹ is the initial value G_(I) ⁰ of the interference signal G_(I), where G_(I) ⁰=[g₁ ¹ g₂ ¹ g₃ ¹ . . . g_(I) ¹].

S107: Repeat steps S103 to S106, until the value X_(S) ¹ of the user signal X_(S) and the value G_(I) ¹ of the interference signal G_(I) that are obtained through an m^(th) iteration meet a preset condition, determine that an approximate value of the user signal X_(S) is X_(S) ¹, and determine that an approximate value of the interference signal G_(I) is G_(I) ¹.

Specifically, the base station performs, by using the foregoing method, iterative processing on the determined value X_(S) ¹ of the user signal X_(S) and the determined value G_(I) ¹ of the interference signal G_(I), until a difference between an estimated error ε^(m-1) obtained through calculation of a formula (4) by using the value X_(S) ¹ of the user signal X_(S) and the value G_(I) ¹ of the interference signal G_(I) that are obtained through an (m-1)^(th) iteration and an estimated error ε^(m) obtained through calculation of the formula (4) by using the value X_(S) ¹ of the user signal X_(S) and the value G_(I) ¹ of the interference signal G_(I) that are obtained through the m^(th) iteration is greater than a preset threshold ε:

ε^(i)=∥Φ_(S) X _(S) ¹+Φ_(I) G _(I) ¹ −y∥₂   (4)

where i is a natural number, and ε^(i) represents an estimated error obtained through an i^(th) iteration; or the base station collects statistics on a quantity of iteration times, until the quantity of iteration times is greater than a preset threshold m.

According to the signal processing method provided in this embodiment, by using respective sparse characteristics of a user signal and an interference signal, the user signal and the interference signal are iteratively recovered, that is, the interference signal is first fastened, and the user signal is recovered; and then the user signal is fastened, and the interference signal is recovered. In this way, an iteration cycle is repeated until the recovered user signal and the recovered interference signal meet a preset condition, thereby effectively improving interference cancellation performance, achieving relatively good interference suppression.

A specific embodiment is used to describe the present invention in detail with reference to a specific application scenario below.

First, the application scenario in this embodiment of the present invention and a process of pre-processing a user signal or an interference signal before the base station receives the signal y are described.

Application Scenario 1

FIG. 2 is a schematic diagram of an application scenario 1 according to the present invention, and an uplink of a cell in a cellular network is considered. It is assumed that K user equipments in the cell communicate with a base station using an OFDM modulation scheme, and a quantity of available subcarriers in a system is N. When user equipment communicates with the base station, the user equipment generally suffers from interference from another system, and for ease of descriptions herein, the another system is referred to as an interference system. Therefore, signals received on a base station side not only include a user signal sent by the user equipment, but also include an interference signal from the another system.

For example, in any timeslot, there are S user equipments that send user signals X_(S) to the base station, and the base station also receives I interference signals G_(I) from the interference system, so that a signal y received by the base station may be represented as:

$y = {{\sum\limits_{k = 1}^{K}{H_{k}{X_{S}(k)}}} + {\sum\limits_{i = 1}^{I}{G_{I}(i)}}}$

where H_(k) is a frequency-domain response of a channel passed through by a signal sent by k^(th) user equipment, X_(S)(k) is a user signal sent by the k^(th) user equipment, and G_(I)(i) is an i^(th) interference signal.

Further, let Φ_(S)=[H₁ H₂ . . . H_(K)], and Φ_(I)=E_(I), where E_(I) is a unit vector,

${X_{S} = \begin{bmatrix} {X_{S}(1)} \\ {X_{S}(2)} \\ \ldots \\ {X_{S}(K)} \end{bmatrix}},{and}$ ${G_{I} = \begin{bmatrix} {G_{I}(1)} \\ {G_{I}(2)} \\ \ldots \\ {G_{I}(I)} \end{bmatrix}},$

so that a formula is obtained:

y=Φ _(S) X _(S)+Φ_(I) G _(I)   (5).

In addition, to fully use a frequency-domain diversity gain and improve performance of interference cancellation and signal recovery, on a user side, a user signal sent by each user equipment is sent after being multiplied by a precoding vector. FIG. 3 is a schematic diagram of precoding processing in the application scenario 1 according to the present invention. Assuming that a precoding vector of a k^(th) user is f_(k), the signal y received by the base station may be represented as:

$y = {{\sum\limits_{k = 1}^{K}{H_{k}f_{k}{X_{S}(k)}}} + {\sum\limits_{i = 1}^{I}{{G_{I}(i)}.}}}$

Further, let Φ_(S)=[H₁f₁ H₂f₂ . . . H_(K)f_(K)], and Φ_(I)=I_(I), where I_(I) is a unit vector,

${X_{S} = \begin{bmatrix} {X_{S}(1)} \\ {X_{S}(2)} \\ \ldots \\ {X_{S}(K)} \end{bmatrix}},{and}$ ${G_{I} = \begin{bmatrix} {G_{I}(1)} \\ {G_{I}(2)} \\ \ldots \\ {G_{I}(I)} \end{bmatrix}},$

so that a formula is obtained:

y=Φ _(S) X _(S)+Φ_(I) G _(I)   (5).

It should be noted that, a precoding vector of a subcarrier of each user meets ∥f_(k)∥=1, and may be randomly generated. For example, a unitary matrix is first randomly generated, and then K columns in the unitary matrix are selected as precoding vectors of K users.

Application Scenario 2

FIG. 4 is a schematic diagram of an application scenario 2 according to the present invention, and an uplink in a cellular network is considered. There are multiple cells in a system, where multiple user equipments in each cell communicate with a base station using an OFDM modulation scheme, and a quantity of available subcarriers in the system is N. In a process in which user equipment in each cell communicates with the base station, the user equipment suffers from interference from user equipment of a neighboring cell, and for ease of descriptions herein, the user equipment of the neighboring cell is referred to as interference user equipment. Therefore, signals received on a base station side not only include a user signal sent by the user equipment, but also include an interference signal sent by the interference user equipment.

It should be noted that, in the following descriptions, only interference cancellation and user signal recovery of a reference base station, for example, a base station 0 in FIG. 3, are considered.

For example, in any timeslot, there are S user equipments that send user signals X_(S) to a base station, and the base station also receives interference signals G_(I) sent from I interference user equipments of a neighboring cell, so that the signal y received by the base station may be represented as:

$y = {{\sum\limits_{k = 1}^{K}{H_{k}{X_{S}(k)}}} + {\sum\limits_{i = 1}^{I}{{\overset{\_}{H}}_{i}{{\overset{\_}{G}}_{I}(l)}}}}$

where H_(k) is a frequency-domain response of a channel passed through by a signal sent by k^(th) user equipment, X_(S)(k) is the signal sent by the k^(th) user equipment, H _(i) is a frequency-domain response of a channel passed through by an interference signal from i^(th) user equipment, and G _(I)(i) is the interference signal from the i^(th) user equipment.

Further, let Φ_(S)=[H₁ H₂ . . . H_(K)], Φ_(I)=[H ₁ H ₂ . . . H _(I)],

${X_{S} = \begin{bmatrix} {X_{S}(1)} \\ {X_{S}(2)} \\ \ldots \\ {X_{S}(K)} \end{bmatrix}},{and}$ ${G_{I} = \begin{bmatrix} {{\overset{\_}{G}}_{I}(1)} \\ {{\overset{\_}{G}}_{I}(2)} \\ \ldots \\ {{\overset{\_}{G}}_{I}(K)} \end{bmatrix}},$

so that a formula is obtained:

y=Φ _(S) X _(S)+Φ_(I) G _(I)   (5).

In addition, similar to the application scenario 1, to fully use a frequency-domain diversity gain and improve performance of interference cancellation and signal recovery, on a user side, a user signal sent by each user equipment is sent after being multiplied by a precoding vector. Assuming that a precoding vector of a k^(th) user is f_(k), the signal y received by the base station may be represented as:

$y = {{\sum\limits_{k = 1}^{K}{H_{k}f_{k}{X_{S}(k)}}} + {\sum\limits_{i = 1}^{I}{{\overset{\_}{H}}_{i}{\overset{\_}{f}}_{i}{{{\overset{\_}{G}}_{I}(i)}.}}}}$

Further, let Φ_(S)=[H₁f₁ H₂f₂ . . . H_(K)f_(K)], Φ_(I)=[H ₁ f ₁ H ₂ f ₂ . . . H _(I) f _(I)],

${X_{S} = \begin{bmatrix} {X_{S}(1)} \\ {X_{S}(2)} \\ \ldots \\ {X_{S}(K)} \end{bmatrix}},{and}$ ${G_{I} = \begin{bmatrix} {{\overset{\_}{G}}_{I}(1)} \\ {{\overset{\_}{G}}_{I}(2)} \\ \ldots \\ {{\overset{\_}{G}}_{I}(K)} \end{bmatrix}},$

so that a formula is obtained:

y=Φ _(S) X _(S)+Φ_(I) G _(I)   (5).

According to the foregoing formula (1), after the base station receives the signal y, recovery of the user signal X_(S) and cancellation of the interference signal X_(I) may be transformed, according to a compressive sensing theory, into calculating approximate solutions of the user signal X_(S) and the interference signal G_(I).

It should be noted that, both Embodiment 2 and Embodiment 3 described below are applicable to the foregoing application scenarios.

Embodiment 2

FIG. 5 is a schematic flowchart of Embodiment 2 of a signal processing method according to the present invention, and the signal processing method specifically includes the following steps:

S201: A base station receives a signal y, where the signal y includes a user signal X_(S) and an interference signal G_(I).

S202: The base station determines an estimated error ε⁰ by using a formula (6):

ε⁰ =∥Φ_(S) X _(S) ⁰+Φ_(I) G _(I) ⁰ −y∥₂   (6)

where ∥ ∥₂ denotes the l₂-norm; Φ_(S)=H or Φ_(S)=H_(S)f_(S), where H_(s) is a frequency-domain response of a channel passed through by the user signal, and f_(s) is a precoding vector used by user equipment UE; Φ_(I)=E_(I) or Φ_(I)=H_(I)f_(I), where E_(I) is a unit vector, H_(I) is a frequency-domain response of a channel passed through by the interference signal, and f_(I) is a precoding vector used by interference user equipment UE; X_(S) ⁰ is an initial value of the user signal X_(S); and G_(I) ⁰ is an initial value of the interference signal G_(I).

Specifically, the base station may preset the initial value X_(S) ⁰ of the user signal X_(S) and the initial value G_(I) ⁰ of the interference signal G_(I), or may obtain through calculation by using the following method: according to the formula (5), further let Φ=[Φ_(S) Φ_(I)], and

${Z = \begin{bmatrix} X_{S} \\ G_{I} \end{bmatrix}},$

and the signal y received by the base station may be further represented as:

y=ΦZ   (7).

The base station transforms, by using a compressive sensing theory, a process of solving the formula (7) into solving the following formula (8) to obtain an optimal solution:

$\begin{matrix} {{\min\limits_{z}{{{\Phi \; Z} - y}}_{2\;}},{{s.t.{Z}_{0}} = {S + I}}} & (8) \end{matrix}$

where S+I is a quantity of nonzero elements of a vector Z, the base station may solve, by using an existing compressive sensing algorithm, for example, a BP algorithm, an MP algorithm, an OMP algorithm, an ROMP algorithm, a CoSaMP algorithm, and an SP algorithm, the formula (8) to obtain Z. It may be known from

$Z = \begin{bmatrix} X_{S} \\ G_{I} \end{bmatrix}$

that, Z includes an optimal solution of the user signal X_(S) and an optimal solution of the interference signal G_(I), and the optimal solution X_(S) ⁰ of the user signal X_(S) and the optimal solution G_(I) ⁰ of the interference signal G_(I) may be obtained by separating the optimal solution of Z, so that it may be determined that X_(S) ⁰ is the initial value of the user signal X_(S), and it may be determined that G_(I) ⁰ is the initial value of the interference signal G_(I). After the initial value of the user signal X_(S) and the initial value of the interference signal G_(I) are determined, the formula (6) is used to calculate an initial estimated error ε⁰.

S203: The base station solves, by using a compressive sensing signal recovery algorithm, a formula (9) to obtain Z_(S), where Z_(S) is an optimal solution X_(S) ¹ of the user signal X_(S):

$\begin{matrix} {{\min\limits_{z_{S}}{{{\Phi_{S}Z_{S}} - y_{S}}}_{2\;}},{{s.t.{Z_{S}}_{0}} = S}} & (9) \end{matrix}$

where min∥ ∥denotes a minimum value; ∥ ∥₀ denotes the l₀-norm; y_(S)=y−Φ_(I)G_(I) ⁰; and S is a quantity of nonzero elements included in the user signal.

Specifically, after the initial value X_(S) ⁰ of the user signal X_(S) and the initial value G_(I) ⁰ of the interference signal G_(I) are determined, the initial value G_(I) ⁰ of the interference signal G_(I) is set to a fixed known value, to continue to calculate the optimal solution of the user signal X_(S). The initial value G_(I) ⁰ of the interference signal G_(I) is substituted into the foregoing formula (5), so as to obtain y=Φ_(S)X_(S)+Φ_(I)G_(I) ⁰, and further obtain Φ_(S)X_(S)=y−Φ_(I)G_(I) ⁰, let y_(S)=y−Φ_(I)G_(I) ⁰, and a formula (10) may be obtained:

y_(S)=Φ_(S)X_(S)   (10).

The compressive sensing theory is used to transform a process of solving the formula (10) into solving the foregoing formula (9) to obtain an optimal solution, and the base station solves, by using the existing compressive sensing signal recovery algorithm, the formula (9) to obtain Z_(S), where Z_(S) is an optimal solution X_(S) ¹ of the user signal X_(S). It should be noted that, a same compressive sensing signal recovery algorithm is used in steps of this embodiment to solve the formula (9) to obtain Z_(S), where Z_(S) is the optimal solution X_(S) ¹ of the user signal X_(S).

S204: The base station solves, by using the compressive sensing signal recovery algorithm, a formula (11) to obtain Z_(I), where Z_(I) is an optimal solution G_(I) ¹ of the interference signal G_(I):

$\begin{matrix} {{\min\limits_{z_{I}}{{{\Phi_{I}Z_{I}} - y_{I}}}_{2\;}},{{s.t.{Z_{I}}_{0}} = I}} & (11) \end{matrix}$

where y_(I)=y−Φ_(I)X_(S) ¹; and I is a quantity of nonzero elements included in the interference signal.

Specifically, similar to step S203, after the optimal solution X_(S) ¹ of the user signal X_(S) is determined, the optimal solution X_(S) ¹ of the user signal X_(S) is set to a fixed known value, to continue to calculate the optimal solution of the interference signal G_(I). The optimal solution X_(S) ¹ of the user signal X_(S) is substituted into the foregoing formula (1), so as to obtain y+Φ_(S)X_(S) ¹+Φ_(I)G_(I), and further obtain Φ_(I)G_(I)=y−Φ_(S)X_(S) ¹, let y−Φ_(S)X_(S) ¹, and a formula (12) may be obtained:

y_(I)=Φ_(I)G_(I)   (12).

The compressive sensing theory is used to transform a process of solving the formula (12) into solving the foregoing formula (11) to obtain an optimal solution, and the base station solves, by using the existing compressive sensing signal recovery algorithm, the formula (12) to obtain Z_(I), where Z_(I) is the optimal solution G_(I) ¹ of the user signal G_(I).

S205: The base station determines an estimated error ε¹ by using a formula (13) to:

ε¹=∥Φ_(S) X _(S) ¹+Φ_(I) G _(I) ¹ −y∥₂   (13).

S206: If a difference between ε⁰ and ε¹ is greater than a preset threshold ε, the base station continues to repeat steps S203 and S204, and solves, by using the compressive sensing signal recovery algorithm, the formula (9) to obtain Z_(S), where Z_(S) is an optimal solution X_(S) ² of the user signal X_(S), where y_(S)=y−Φ_(I)G_(I) ^(I); the base station solves, by using the compressive sensing signal recovery algorithm, the formula (11) to obtain Z_(I), where Z_(I) is an optimal solution G_(I) ² of the interference signal G_(I), where y_(I)=y−Φ_(I)X_(S) ²; and the base station determines an estimated error ε² by using the formula (10):

ε²=∥Φ_(S) X _(S) ²+Φ_(I) G _(I) ² −y∥₂   (14).

If a difference between ε² and ε¹ is greater than the preset threshold ε, an iteration process in steps S203 and S204 is repeated, until a difference between ε^(m-1) and ε^(m) that are obtained by using an optimal solution X_(S) ^(m) of a sent signal X_(S) and an optimal solution G_(I) ^(m) of the interference signal G_(I) that are obtained through an m^(th) calculation is less than or equal to the preset threshold ε, where X_(S) ^(m) is an estimated value of the user signal X_(S), and G_(I) ^(m) is an estimated value of the interference signal G_(I).

In this embodiment, by using respective sparse characteristics of a user signal and an interference signal, the user signal and the interference signal are iteratively recovered by using a compressive sensing signal recovery algorithm, that is, the interference signal is first fastened, and the user signal is recovered by using the compressive sensing signal recovery algorithm; and then the user signal is fastened, and the interference signal is recovered by using the compressive sensing signal recovery algorithm. In this way, an iteration cycle is repeated until the recovered user signal and the recovered interference signal meet a preset error, thereby effectively improving interference cancellation performance, achieving relatively good interference suppression.

In addition, a joint interference signal cancellation and compressive sensing signal recovery algorithm based on sparseness maximum likelihood estimation is used in this embodiment. In a case without knowing a user signal and an interference signal, remaining signal latitude that is not occupied by the user signal and the interference signal can be fully used, to improve interference signal cancellation performance, and a signal-to-noise ratio and a diversity gain that are of a receive end. Because joint interference signal cancellation and user signal recovery is used, and the respective sparse characteristics of the user signal and the interference signal are used, interference cancellation performance in this embodiment is greatly better than that in an existing interference cancellation method, and none of the following limitation conditions is needed in this embodiment: 1. a system has extra latitude; 2. interference is a sparse frequency domain; 3. support of interference and support of a wanted signal are not overlapped; or 4. the support of the interference is known. Therefore, an application scope of this embodiment is wider.

Embodiment 3

FIG. 6 is a schematic flowchart of Embodiment 3 of a signal processing method according to the present invention, and a difference between this embodiment and Embodiment 2 is that in this embodiment, a quantity of iteration times for calculating a user signal and an interference signal is controlled by using a quantity of cycles. The method specifically includes the following steps:

S301: A base station receives a signal y, where the signal y includes a user signal X_(S) and an interference signal G_(I).

S302: The base station initializes a quantity n of cycles to 0, and solves, by using a compressive sensing signal recovery algorithm, a formula (15) to obtain Z_(S), where Z_(S) is an optimal solution X_(S) ¹ of the user signal X_(S):

$\begin{matrix} {{\min\limits_{z_{S}}{{{\Phi_{S}Z_{S}} - y_{S}}}_{2\;}},{{s.t.{Z_{S}}_{0}} = S}} & (15) \end{matrix}$

where min∥ ∥ denotes a minimum value; ∥ ∥₀ denotes the l₀-norm; ∥ ∥₂ denotes the l₂-norm; Φ_(S)=H or Φ_(S)=H_(S)f_(S), where H_(s) is a frequency-domain response of a channel passed through by the user signal, and f_(s) is a precoding vector used by user equipment UE; y_(S)=y−Φ_(I)G_(I) ⁰; S is a quantity of nonzero elements included in the user signal; and G_(I) ⁰ is an initial value of the interference signal G_(I).

Specifically, similar to step S202 in Embodiment 1, the base station may preset an initial value G_(I) ⁰ of the interference signal G_(I), or may obtain through calculation by using the following method: according to the formula (5), further let Φ=[Φ_(S) Φ_(I)], and

${Z = \begin{bmatrix} X_{S} \\ G_{I} \end{bmatrix}},$

and the signal y received by the base station may be further represented as:

y=ΦZ   (7).

The base station transforms, by using a compressive sensing theory, a process of solving the formula (7) into solving the following formula (8) to obtain an optimal solution:

$\begin{matrix} {{\min\limits_{Z}{{{\Phi \; Z} - y}}_{2}},{{s.t.{Z}_{0}} = {S + I}}} & (8) \end{matrix}$

where S+I is a quantity of nonzero elements of a vector Z, the base station may solve, by using an existing compressive sensing algorithm, for example, a BP algorithm, an MP algorithm, an OMP algorithm, an ROMP algorithm, a CoSaMP algorithm, and an SP algorithm, the formula (8) to obtain Z. It may be known from

$Z = \begin{bmatrix} X_{S} \\ G_{I} \end{bmatrix}$

that, Z includes an optimal solution of the user signal X_(S) and an optimal solution of the interference signal G_(I), and the optimal solution G_(I) ⁰ of the interference signal G_(I) may be obtained by separating an optimal solution of Z, so that it may be determined that G_(I) ⁰ is the initial value of the interference signal G_(I). After the initial value G_(I) ⁰ of the interference signal G_(I) is determined, the initial value G_(I) ⁰ of the interference signal G_(I) is set to a fixed known value, to calculate the optimal solution of the user signal X_(S). The initial value G_(I) ⁰ of the interference signal G_(I) is substituted into the foregoing formula (5), so as to obtain y=Φ_(S)X_(S)+Φ_(I)G_(I) ⁰, and further obtain Φ_(S)X_(S)=y−Φ_(I)G_(I) ⁰. Let y_(S)=y−Φ_(I)G_(I) ⁰, and the formula (10) may be obtained. A compressive sensing theory is used to transform a process of solving the formula (10) into solving the foregoing formula (9) to obtain an optimal solution. The base station solves, by using the existing compressive sensing signal recovery algorithm, the formula (9) to obtain Z_(S), where Z_(S) is an optimal solution X_(S) ¹ of the user signal X_(S).

S303: The base station solves, by using the compressive sensing signal recovery algorithm, a formula (16) to obtain Z_(I), where Z_(I) is an optimal solution G_(I) ¹ of the interference signal G_(I):

$\begin{matrix} {{\min\limits_{Z_{I}}{{{\Phi_{I}Z_{I}} - y_{I}}}_{2}},{{s.t.{Z_{I}}_{0}} = I}} & (16) \end{matrix}$

where Φ_(I)=E or Φ_(I)=H_(I)f_(I), where E is a unit vector, H_(I) is a frequency-domain response of a channel passed through by the interference signal, and f_(I) is a precoding vector used by interference user equipment UE; y_(I)=y−Φ_(I)X_(S) ¹; and I is a quantity of nonzero elements included in the interference signal.

For details, refer to step S203 in Embodiment 1, and details are not described herein again.

S304: The base station updates the quantity n of cycles to 1, and solves, by using the compressive sensing signal recovery algorithm, the formula (15) to obtain Z_(S), where Z_(S) is an optimal solution X_(S) ² of the user signal X_(S), where y_(S)=y−Φ_(I)G_(I) ¹.

S305: The base station solves, by using the compressive sensing signal recovery algorithm, the formula (16) to obtain Z_(I), where Z_(I) is an optimal solution G_(I) ² of the interference signal G_(I), where y_(I)=y−Φ_(I)X_(S) ².

S306: If n is less than or equal to a preset threshold m, the base station repeats an iteration process of steps S304 and S305, until the quantity n of cycles is greater than the preset threshold m, where an optimal solution X_(S) ^(m) of a sent signal X_(S) obtained through an m^(th) calculation is an estimated value of the user signal X_(S), and an optimal solution G_(I) ^(m) of the interference signal G_(I) obtained through the m^(th) calculation is an estimated value of the interference signal G_(I).

In this embodiment, by using respective sparse characteristics of a user signal and an interference signal, the user signal and the interference signal are iteratively recovered by using a compressive sensing signal recovery algorithm, that is, the interference signal is first fastened, and the user signal is recovered by using the compressive sensing signal recovery algorithm; and then the user signal is fastened, and the interference signal is recovered by using the compressive sensing signal recovery algorithm. In this way, an iteration cycle is repeated until the recovered user signal and the recovered interference signal meet a quantity of cycles, thereby effectively improving interference cancellation performance, achieving relatively good interference suppression.

Embodiment 3

FIG. 7 is a schematic structural diagram of an embodiment of a base station according to the present invention, and the base station specifically includes:

a receiving module 41, configured to receive a signal y, where the signal y includes a user signal X_(S) and an interference signal G_(I), where: X_(S) is a vector including S user signals, and X_(S)=[x₁ x₂ x₃ . . . x_(s)], where x₁ represents a first user signal, x₂ represents a second user signal, x₃ represents a third user signal, and so on, and x_(s) represents an s^(th) user signal; and G_(I)is a vector including I interference signals, and G_(I)=[g₁ g₂ g₃ . . . g_(I)], where g₁ represents a first interference signal, g₂ represents a second interference signal, g₃ represents a third interference signal, and so on, and g_(I) represents an interference signal; and

a determining module 42, configured to determine an initial value G_(I) ⁰ of the interference signal G_(I), where G_(I) ⁰=[g₁ ⁰ g₂ ⁰ g₃ ⁰ . . . g_(I) ⁰]; determine a value X_(S) ¹ of the user signal X_(S) according to the initial value G_(I) ⁰ and the signal y, where X_(S) ¹=[x₁ ¹ x₂ ¹ x₃ ¹ . . . x_(s) ¹]; determine that the value X_(S) ¹ is an initial value X_(S) ⁰ of the user signal X_(S), where X_(S) ⁰=[x₁ ¹ x₂ ¹ x₃ ¹ . . . x_(s) ¹]; determine a value G_(I) ¹ of the interference signal G_(I) according to the initial value X_(S) ⁰ and the signal y, where G_(I) ¹=[g₁ ¹ g₂ ¹ g₃ ¹ . . . g_(I) ¹]; determine that the value G_(I) ¹ is the initial value G_(I) ⁰ of the interference signal G_(I), where G_(I) ⁰=[g_(I) ¹ g₂ ¹ g₃ ¹ . . . g_(I) ¹]; and perform iterative processing on the value X_(S) ¹ of the user signal X_(S) and the value G_(I) ¹ of the interference signal G_(I), until the value X_(S) ¹ of the user signal X_(S) and the value G_(I) ¹ of the interference signal G_(I) that are obtained through an m^(th) iteration meet a preset condition, determine that an approximate value of the user signal X_(S) is X_(S) ¹, and determine that an approximate value of the interference signal G_(I) is G_(I) ¹.

Further, the determining module 42 is specifically configured to:

solve a formula (1) to obtain Z, where Z includes an optimal solution G_(I) ^(opt) of the interference signal G_(I):

$\begin{matrix} {{\min\limits_{Z}{{{\Phi \; Z} - y}}_{2}},{{s.t.{Z}_{0}} = {S + I}}} & (1) \end{matrix}$

where ∥ ∥₂ denotes the l₂-norm; ∥ ∥₀ denotes the l₀-norm; Φ=[Φ_(S) Φ_(I)], and

${Z = \begin{bmatrix} X_{S} \\ G_{I} \end{bmatrix}};$

Φ_(S)=H or Φ_(S)=H_(S)f_(S), where H_(s) is a frequency-domain response of a channel passed through by the user signal, and f_(s) is a precoding vector used by user equipment UE; and Φ_(I)=E_(I) or Φ_(I)=H_(I)f_(I), where E_(I) is a unit vector, H_(I) is a frequency-domain response of a channel passed through by the interference signal, and f_(I) is a precoding vector used by interference user equipment UE; and determine that the optimal solution G_(I) ^(opt) is the initial value G_(I) ⁰ of the interference signal G_(I).

Further, the determining module 42 is specifically configured to:

solve a formula (2) to obtain an optimal solution Z_(S), and determine that Z_(S) is the value X_(S) ¹ of the user signal X_(S):

$\begin{matrix} {{\min\limits_{Z_{S}}{{{\Phi_{S}Z_{S}} - y_{S}}}_{2}},{{s.t.{Z_{S}}_{0}} = S}} & (2) \end{matrix}$

where y_(S)=y−Φ_(I)G_(I) ⁰; and solve a formula (3) to obtain an optimal solution Z_(I), and determine that Z_(I) is the value G_(I) ¹ of the interference signal G_(I):

$\begin{matrix} {{\min\limits_{Z_{I}}{{{\Phi_{I}Z_{I}} - y_{I}}}_{2}},{{s.t.{Z_{I}}_{0}} = I}} & (3) \end{matrix}$

where y_(I)=y−Φ_(S)X_(S) ⁰.

Further, the base station further includes: a compressive sensing module 43, configured to solve, by using a compressive sensing signal recovery algorithm, a formula to obtain an optimal solution.

Further, the determining module 42 is specifically configured to:

perform iterative processing on the value X_(S) ¹ of the user signal X_(S) and the value G_(I) ¹ of the interference signal G_(I), until a difference between an estimated error ε^(m-1) obtained through calculation of a formula (4) by using the value X_(S) ¹ of the user signal X_(S) and the value G_(I) ¹ of the interference signal G_(I) that are obtained through an (m-1)^(th) iteration and an estimated error ε^(m) obtained through calculation of the formula (4) by using the value X_(S) ¹ of the user signal X_(S) and the value G_(I) ¹ of the interference signal G_(I) that are obtained through the m^(th) iteration is greater than a preset threshold ε:

ε^(i)=∥Φ_(S) X _(S) ¹+Φ_(I) G _(I) ¹ −y∥₂   (4)

where i is a natural number, and ε^(i) represents an estimated error obtained through an i^(th) iteration; or collect statistics on a quantity of iteration times, until the quantity of iteration times is greater than a preset threshold m.

The base station in this embodiment is configured to execute the solutions of the method embodiments in Embodiment 1, Embodiment 2, and Embodiment 3. The implementation principle and the technical effect of this embodiment are similar to those of the method embodiments, and are not described in detail herein.

Persons of ordinary skill in the art may understand that all or some of the steps of the method embodiments may be implemented by a program instructing relevant hardware. The program may be stored in a computer-readable storage medium. When the program runs, the steps of the method embodiments are performed. The foregoing storage medium includes: any medium that can store program code, such as a ROM, a RAM, a magnetic disk, or an optical disc.

Finally, it should be noted that the foregoing embodiments are merely intended for describing the technical solutions of the present invention, but not for limiting the present invention. Although the present invention is described in detail with reference to the foregoing embodiments, persons of ordinary skill in the art should understand that they may still make modifications to the technical solutions described in the foregoing embodiments or make equivalent replacements to some or all technical features thereof, without departing from the scope of the technical solutions of the embodiments of the present invention. 

What is claimed is:
 1. A signal processing method, comprising: receiving, by a base station, a signal y, wherein the signal y comprises a user signal X_(S) and an interference signal G_(I), wherein: X_(S) is a vector comprising S user signals, and X_(S)=[x₁ x₂ x₃ . . . x_(s)], wherein x₁ represents a first user signal, x₂ represents a second user signal, x₃ represents a third user signal, and so on, and x_(s) represents an s^(th) user signal; and G_(I) is a vector comprising I interference signals, and G_(I)=[g₁ g₂ g₃ . . . g_(I)], wherein g₁ represents a first interference signal, g₂ represents a second interference signal, g₃ represents a third interference signal, and so on, and g_(I) represents an I^(th) interference signal; determining, by the base station, an initial value G_(I) ⁰ of the interference signal G_(I), wherein G_(I) ⁰=[g₁ ⁰ g₂ ⁰ g₃ ⁰ . . . g_(I) ⁰]; determining, by the base station, a value X_(S) ¹ of the user signal X_(S) according to the initial value G_(I) ⁰ and the signal y, wherein X_(S) ¹=[x₁ ¹ x₂ ¹ x₃ ¹ . . . x_(s) ¹]; determining, by the base station, that the value X_(S) ¹ is an initial value X_(S) ⁰ of the user signal X_(S), wherein X_(S) ⁰=[x₁ ¹ x₂ ¹ x₃ ¹ . . . x_(s) ¹]; determining, by the base station, a value G_(I) ¹ of the interference signal G_(I) according to the initial value X_(S) ⁰ and the signal y, wherein G_(I) ¹=[g₁ ¹ g₂ ¹ g₃ ¹ . . . g_(I) ¹]; determining, by the base station, that the value G_(I) ¹ is the initial value G_(I) ⁰ of the interference signal G_(I), wherein G_(I) ⁰=[g₁ ¹ g₂ ¹ g₃ ¹ . . . g_(I) ¹]; and performing, by using the foregoing method, iterative processing on the determined value X_(S) ¹ of the user signal X_(S) and the determined value G_(I) ¹ of the interference signal G_(I), until the value X_(S) ¹ of the user signal X_(S) and the value G_(I) ¹ of the interference signal G_(I) that are obtained through an m^(th) iteration meet a preset condition, determining that an approximate value of the user signal X_(S) is X_(S) ¹, and determining that an approximate value of the interference signal G_(I) is G_(I) ¹.
 2. The method according to claim 1, wherein the determining, by the base station, an initial value G_(I) ⁰ of the interference signal G_(I) comprises: solving, by the base station, a formula (1) to obtain Z, wherein Z comprises an optimal solution G_(I) ^(opt) of the interference signal G_(I): $\begin{matrix} {{\min\limits_{Z}{{{\Phi \; Z} - y}}_{2}},{{s.t.{Z}_{0}} = {S + I}}} & (1) \end{matrix}$ wherein ∥ ∥₂ denotes the l₂-norm; ∥ ∥₀ denotes the l₀-norm; Φ=[Φ_(S) Φ_(I)], and ${Z = \begin{bmatrix} X_{S} \\ G_{I} \end{bmatrix}};$ Φ_(S)=H or Φ_(S)=H_(S)f_(S), wherein H_(s) is a frequency-domain response of a channel passed through by the user signal, and f_(s) is a precoding vector used by user equipment UE; and Φ_(I)=E_(I) or Φ_(I)=H_(I)f_(I), wherein E_(I) is a unit vector, H_(I) is a frequency-domain response of a channel passed through by the interference signal, and f_(I) is a precoding vector used by interference user equipment UE; and determining, by the base station, that the optimal solution G_(I) ^(opt) is the initial value G_(I) ⁰ of the interference signal G_(I).
 3. The method according to claim 1, wherein the determining, by the base station, a value X_(S) ¹ of the user signal X_(S) according to the initial value G_(I) ⁰ and the signal y comprises: solving, by the base station, a formula (2) to obtain an optimal solution Z_(S), and determining that Z_(S) is the value X_(S) ¹ of the user signal X_(S): $\begin{matrix} {{\min\limits_{Z_{S}}{{{\Phi_{S}Z_{S}} - y_{S}}}_{2}},{{s.t.{Z_{S}}_{0}} = S}} & (2) \end{matrix}$ wherein y_(S)=y−Φ_(I)G_(I) ⁰; and correspondingly, the determining, by the base station, a value G_(I) ¹ of the interference signal G_(I) according to the initial value X_(S) ⁰ and the signal y comprises: solving, by the base station, a formula (3) to obtain an optimal solution Z_(I), and determining that Z_(I) is the value G_(I) ¹ of the interference signal G_(I): $\begin{matrix} {{\min\limits_{Z_{I}}{{{\Phi_{I}Z_{I}} - y_{I}}}_{2}},{{s.t.{Z_{I}}_{0}} = I}} & (3) \end{matrix}$ wherein y_(I)=y−Φ_(S)X_(S) ⁰.
 4. The method according to claim 2, wherein the base station solves a formula by using a compressive sensing signal recovery algorithm to obtain an optimal solution.
 5. The method according to claim 1, wherein the performing, by the base station by using the foregoing method, iterative processing on the determined value X_(S) ¹ of the user signal X_(S) and the determined value G_(I) ¹ of the interference signal G_(I), until the value X_(S) ¹ of the user signal X_(S) and the value G_(I) ¹ of the interference signal G_(I) that are obtained through an m^(th) iteration meet a preset condition comprises. performing, by the base station by using the foregoing method, iterative processing on the determined value X_(S) ¹ of the user signal X_(S) and the determined value G_(I) ¹ of the interference signal G_(I), until a difference between an estimated error ε^(m-1) obtained through calculation of a formula (4) by using the value X_(S) ¹ of the user signal X_(S) and the value G_(I) ¹ of the interference signal G_(I) that are obtained through an (m-1)^(th) iteration and an estimated error ε^(m) obtained through calculation of the formula (4) by using the value X_(S) ¹ of the user signal X_(S) and the value G_(I) ¹ of the interference signal G_(I) that are obtained through the m^(th) iteration is greater than a preset threshold ε: ε^(i)=∥Φ_(S) X _(S) ¹+Φ_(I) G _(I) ¹ −y∥₂   (4) wherein i is a natural number, and ε^(i) represents an estimated error obtained through an i^(th) iteration; or collecting, by the base station, statistics on a quantity of iteration times, until the quantity of iteration times is greater than a preset threshold m.
 6. A base station, comprising: a receiving module, configured to receive a signal y, wherein the signal y comprises a user signal X_(S) and an interference signal G_(I), wherein: X_(S) is a vector comprising S user signals, and X_(S) =[x₁ x₂ x₃ . . . x_(s)], wherein x₁ represents a first user signal, x₂ represents a second user signal, x₃ represents a third user signal, and so on, and x_(s) represents an S^(th) user signal; and G_(I) is a vector comprising I interference signals, and G_(I)=[g₁ g₂ g₃ . . . g_(I)], wherein g₁ represents a first interference signal, g₂ represents a second interference signal, g₃ represents a third interference signal, and so on, and g_(I) represents an I^(th) interference signal; and a determining module, configured to determine an initial value G_(I) ⁰ of the interference signal G_(I), wherein G_(I) ⁰=[g₁ ⁰ g₂ ⁰ g₃ ⁰ . . . g_(I) ⁰]; determine a value X_(S) ¹ of the user signal X_(S) according to the initial value G_(I) ⁰ and the signal y, wherein X_(S) ¹=[x₁ ¹ x₂ ¹ x₃ ¹ . . . x_(s) ¹]; determine that the value X_(S) ¹ is an initial value X_(S) ⁰ of the user signal X_(S), wherein X_(S) ⁰=[x₁ ¹ x₂ ¹ x₃ ¹ . . . x_(s)]; determine a value G_(I) ¹ of the interference signal G_(I) according to the initial value X_(S) ⁰ and the signal y, wherein G_(I) ¹=[g₁ ¹ g₂ ¹ g₃ ¹ . . . g_(I) ¹]; determine that the value G_(I) ¹ is the initial value G_(I) ⁰ of the interference signal G_(I), wherein G_(I) ⁰=[g₁ ¹ g₂ ¹ g₃ ¹ . . . g_(I) ¹]; and perform iterative processing on the value X_(S) ¹ of the user signal X_(S) and the value G_(I) ¹ of the interference signal G_(I), until the value X_(S) ¹ of the user signal X_(S) and the value G_(I) ¹ of the interference signal G_(I) that are obtained through an m^(th) iteration meet a preset condition, determine that an approximate value of the user signal X_(S) is X_(S) ¹, and determine that an approximate value of the interference signal G_(I) is G_(I) ¹.
 7. The base station according to claim 6, wherein the determining module is specifically configured to: solve a formula (1) to obtain Z, wherein Z comprises an optimal solution G_(I) ^(opt) of the interference signal G_(I): $\begin{matrix} {{\min\limits_{Z}{{{\Phi \; Z} - y}}_{2}},{{s.t.{Z}_{0}} = {S + I}}} & (1) \end{matrix}$ wherein ∥ ∥₂ denotes the l₂-norm; ∥ ∥₀ denotes the l₀-norm; Φ=[Φ_(S) Φ_(I)], and ${Z = \begin{bmatrix} X_{S} \\ G_{I} \end{bmatrix}};$ Φ_(S)=H or Φ_(S)=H_(S)f_(S), wherein H_(s) is a frequency-domain response of a channel passed through by the user signal, and f_(s) is a precoding vector used by user equipment UE; and Φ_(I)=E_(I) or Φ_(I)=H_(I) f_(I), wherein E_(I) is a unit vector, H_(I) is a frequency-domain response of a channel passed through by the interference signal, and f_(I) is a precoding vector used by interference user equipment UE; and determine that the optimal solution G_(I) ^(opt) is the initial value G_(I) ⁰ of the interference signal G_(I).
 8. The base station according to claim 6, wherein the determining module is specifically configured to: solve a formula (2) to obtain an optimal solution Z_(S), and determine that Z_(S) is the value X_(S) ¹ of the user signal X_(S): $\begin{matrix} {{\min\limits_{Z_{S}}{{{\Phi_{S}Z_{S}} - y_{S}}}_{2}},{{s.t.{Z_{S}}_{0}} = S}} & (2) \end{matrix}$ wherein y_(S)=y−Φ_(I)G_(I) ⁰; and solve a formula (3) to obtain an optimal solution Z_(I), and determine that Z_(I) is the value G_(I) ¹ of the interference signal G_(I): $\begin{matrix} {{\min\limits_{Z_{I}}{{{\Phi_{I}Z_{I}} - y_{I}}}_{2}},{{s.t.{Z_{I}}_{0}} = I}} & (3) \end{matrix}$ wherein y_(I)=y−Φ_(S)X_(S) ⁰.
 9. The base station according to claim 7, further comprising a compressive sensing module, configured to solve, by using a compressive sensing signal recovery algorithm, a formula to obtain an optimal solution.
 10. The base station according to claim 6, wherein the determining module is specifically configured to: perform iterative processing on the value X_(S) ¹ of the user signal X_(S) and the value G_(I) ¹ of the interference signal G_(I), until a difference between an estimated error ε^(m-1) obtained through calculation of a formula (4) by using the value X_(S) ¹ of the user signal X_(S) and the value G_(I) ¹ of the interference signal G_(I) that are obtained through an (m-1)^(th) iteration and an estimated error ε^(m) obtained through calculation of the formula (4) by using the value X_(S) ¹ of the user signal X_(S) and the value G_(I) ¹ of the interference signal G_(I) that are obtained through the m^(th) iteration is greater than a preset threshold ε: ε^(i)=∥Φ_(S) X _(S) ¹+Φ_(I) G _(I) ¹ −y∥₂   (4) wherein i is a natural number, and ε^(i) represents an estimated error obtained through an i^(th) iteration; or collect statistics on a quantity of iteration times, until the quantity of iteration times is greater than a preset threshold m. 